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Abstract 



(N 

O 
<N 

, This article deals with the existence and the uniqueness of solutions to quadratic and superquadratic Marko- 

vian backward stochastic differential equations (BSDEs for short) with an unbounded terminal condition. Our 
results are deeply linked with a strong a priori estimate on Z that takes advantage of the Markovian framework. 
This estimate allows us to prove the existence of a viscosity solution to a semilinear parabolic partial differen- 
tial equation with nonlinearity having quadratic or superquadratic growth in the gradient of the solution. This 
, estimate also allows us to give explicit convergence rates for time approximation of quadratic or superquadratic 

p I . Markovian BSDEs. 

■S: 

" 1 Introduction 

Since the early nineties and the work of Pardoux and Peng Il25l , there has been an increasing interest for backward 
CSJ ■ stochastic differential equations (BSDEs for short). These equations have a wide range of applications in stochastic 

, control, in finance or in partial differential equation theory. A particular class of BSDE is studied since few years: 

BSDEs with generators of quadratic growth with respect to the variable z (quadratic BSDEs for short). This 
class arises, for example, in the context of utility optimization problems with exponential utility functions, or 
ly-^ ■ alternatively in questions related to risk minimization for the entropic risk measure (see e.g. Il27l fTTl |24| among 

many other references). Many papers deal with existence and uniqueness of solution for such BSDEs. In the first 
one II2TI . Kobylanski obtains an existence and uniqueness result for quadratic BSDEs when the terminal condition 
is bounded. Let us remark that this result has been revisited recently thanks to a fixed point argument by Tevzadze 
in ll28l . Now, it is well known that the boundedness of the terminal condition is a too strong assumption. Indeed, 
when we look to the simple quadratic BSDE 



T \^ |2 



Jt ^ Jt 

we find the explicit solution = log (E [e^|J-t]) and we immediately see that we just need to have an exponential 
moment for ^ to obtain a solution. In fS], Briand and Hu show an existence result for quadratic BSDEs when the 
terminal condition has such an assumption. Let us notice that this result has been recently revisited in fSl by a 
direct forward method that does not use the result of Kobylanski. In this paper, Barrieu and El Karoui obtain 
a monotone stability result for general quadratic semimartingales and then derive an existence result for general 
quadratic BSDEs. For the uniqueness problem, results are more incomplete. In ifTSl . authors show a uniqueness 
result when the generator is convex (or concave) with respect to z and when ^ has an exponential moment which 
is almost the exponential moment needed for the existence result. 

Naturally, we could also wonder what happens when the generator has a superquadratic growth with respect to 
the variable z. Up to our knowledge the case of superquadratic BSDEs is only investigate in the recent paper llT2ll . 
In this article, authors consider superquadratic BSDEs when the terminal condition is bounded and the generator 
is convex in z. Firstly, they show that in a general way the problem is ill-posed: given a superquadratic generator, 
there exists a bounded terminal condition such that the associated BSDE does not admit any bounded solution 
and, on the other hand, if the BSDE admits a bounded solution, there exist infinitely many bounded solutions for 
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this BSDE. In the same paper, authors also show that the problem becomes well-posed in a Markovian framework: 
When the terminal condition and the generator are deterministic functions of a forward SDE, we have an existence 
result. 

The first aim of this paper is to study existence and uniqueness results for quadratic and superquadratic Marko- 
vian BSDEs. More precisely, we consider {X, Y, Z) the solution to the (decoupled) forward backward system 



where / has quadratic or superquadratic growth with respect z, has no convexity assumption, and g is not supposed 
to be bounded. The starting point of our work is a simple result that says; if g and / are Lipschitz functions with 
respect to x, then there exists a unique solution such that Z is bounded, or in other words, Z preserves the regularity 
of the derivatives of g and / with respect to x. Now, the idea is to show that this property stays true when g and / 
are only locally Lipschitz. More precisely, if we assume that 



Thanks to this kind of estimate, it is then possible to show an existence and uniqueness result amongst solutions 
that, roughly speaking, verify such an estimate (see Theorem 12.5b . Contrarily to ITSl ITSll we do not need a 
convexity assumption on / and contrarily to lfT2l . we treat the case of unbounded terminal conditions. On the 
other hand, for the quadratic case we need a framework which is far more restrictive than the general framework 
of lfT3l 121 because we only consider markovian BSDEs with assumptions on the derivatives of g and / instead of 
classical assumptions on the growth of g and /. 

One of the major drawback of results explained before is that we consider only the case of a deterministic 
function a. The second part of our paper gives some partial results when a is random. In this framework we do 
not know if our previous starting point stays true: if g and / are Lipschitz functions with respect to x and if a is 
bounded, does there exist a solution such that Z is bounded ? We are able to show that this is true when T is small 
enough or for all T when we consider a simple example of quadratic BSDE. But the general case stays an open 
question. We also investigate precisely the quadratic case when g and / are bounded with respect to x by deeply 
using bounded mean oscillation martingale (BMO martingale for short) tools. 

Thanks to our existence and uniqueness result we are able to give a probabilistic representation of the following 



Such a probabilistic representation, also called Feynman-Kac representation, is already given in lfT3l when / has 
a quadratic growth and is convex with respect to z. Existence and uniqueness of this PDE has been studied in 
ll9) when / has a quadratic growth with respect to Vucr and in ifTOl for the superquadratic case, but the main part 
of the results needs a convexity assumption on / with respect to z. In this paper, our existence result arises in 
quadratic and superquadratic frameworks. Moreover, we do not need any convexity assumption on / which is 
interesting for applications: For example, when we consider Isaacs equations in differential game theory, / is the 
sum of a convex function and a concave function with respect to z. 

The main goal of this paper is to apply a priori estimates obtained for the process Z to the problem of time 
discretization of quadratic and superquadratic BSDEs. Actually, the design of efficient algorithms which are able 
to solve BSDEs in any reasonable dimension has been intensively studied since the first work of Chevance JS], see 
for instance ||29l [3l lT6l . But in all these works, the driver of the BSDE is a Lipschitz function with respect to z and 
this assumption plays a key role in theirs proofs. In a recent paper, Cheridito and Stadje Q study approximation of 
BSDEs by backward stochastic difference equations which consist in replacing the Brownian motion by a random 
walk. They obtain a convergence result when the driver has a subquadratic growth with respect to z and they give 
an example where there proof does not work when the driver has a quadratic growth. To the best of our knowledge, 
the only works where the time approximation of a quadratic BSDE is studied are the one of Imkeller and dos 
Reis ifTHl and the one of Richou E6i . Let us notice that, when the driver has a specific forrrQ, it is possible to get 

' Roughly speaking, the driver is a sum of a quadratic term z C \ z\^ and a function that has a hnear growth with respect to z. 




|V.g(a;)| + |V,/(.,x,.,.)KC(l + |xr') 



for r sufficiently small, then we are able to show the a priori estimate 



\z\^c{i + \xn. 



PDE: 



dtu{t,x) + £u{t,x) + f{t,x,u{t,x)*Wu{t,x)a{t)) =0, x eR'^,te [0,T], 
u{T, .) = g. 
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around the problem by using an exponential transformation method (see |fT9l ) or by using results on fully coupled 
forward-backward differential equations (see flT]). Papers ifTSl l26l only study the case of a bounded terminal 
condition; The first one investigates the case of Lipschitz terminal conditions whereas the second one studies the 
non-smooth case. To the best of our knowledge, the time approximation of superquadratic BSDEs has not been 
studied yet. In this paper we have obtained two types of results. Firstly we consider the case of a deterministic 
function a. Theorem lS. 71 gives us a speed of convergence very close to the speed of convergence in the classical 
Lipschitz case and this theorem is obtained in a general framework (quadratic and superquadratic BSDEs with 
an unbounded terminal condition). When a is random, we only study quadratic BSDEs with bounded terminal 
conditions. In Theorem lS. 9l we obtain almost the classical speed of convergence but in a restricted framework that 
does not cover some interesting situations: for example we are not able to find a "good" speed of convergence 
when cr and g are Lipschitz functions with respect to x and this question is actually a real challenge. 

The paper is organized as follows. In section 2 we obtain an existence and uniqueness result and an a priori 
estimate on Z for Markovian quadratic and superquadratic BSDEs with unbounded terminal conditions when a 
is a deterministic function. In section 3 we give some extra partial results when cr is random. Section 4 contains 
an application to semilinear parabolic PDEs. The last section is devoted to time approximation of quadratic and 
superquadratic Markovian BSDEs. 



Notations Throughout this paper, {Wt)t^o will denote a d-dimensional Brownian motion, defined on a proba- 
bility space (O, J^, P). For t ^ 0, let J"t denotes the cr-algebra a{Ws;0 ^ s ^ t), augmented with the P-null 
sets of The Euclidean norm on R'' will be denoted by | . | . The operator norm induced by | . | on the space of 
linear operators is also denoted by | . | . The notation Et stands for the conditional expectation given Tf. For p ^ 2, 
TO G N, we denote further 

• iS''(M™), or SP when no confusion is possible, the space of all adapted processes {Yt)t£io.T] with values 
in M™ normed by \\Y\\g^ = E[(supjg[o \Yt\)P]^/P; 5°°(R™), or S°°, the space of bounded measurable 
processes; 

• A^f'(R™), or AiP, the space of all progressively measurable processes (^t)tG[o.T] with values in R™ 
normed by \\Z\\j^, = EK/J' \zf dsr/^]^/r> . 

In the following, we keep the same notation C for all finite, nonnegative constants that appear in our computations. 
In this paper we will consider X the solution to the SDE 

Xt=x+ f bis,Xs)ds+ f a{s)dWs, (1.1) 
Jo Jo 

and (Y, Z) e S'^ x the solution to the Markovian BSDE 

Yt=g{XT)+ [ f{s,Xs,Y,,Zs)ds~ [ Z.dW^. (1.2) 
Jt Jt 

2 A uniqueness and existence result 

For the SDE dl.lb we use standard assumption. 

Assumption (F.l). Let b : [0, T] x R'^ — > R'' and a : [0, T] — > R''^'* be continuous functions and let us assume 
that there exists Kb ^ such that: 

1. yt e [0,T], \b{t,o)\ C, 

2. Vt e [0,T],V(a;,a;') G R'^ x R'^, \b{t,x) ~ b{t,x')\ ^ Kb \x~x'\ . 

Let us assume the following for the generator and the terminal condition of the BSDE il.2l . 

Assumption (B.l). Let / : [0, T] x M'' x R x R^^'' ^ R and g : R'' ^ M be continuous functions and let us 
assume moreover that there exist five constants, 1^1, a^O, /3^0, j^O and Kfy ^ such that: 

1. for each {t,x,y,y',z) G [0, T] x R'' x R x R x M.^""^, 



\f{t,x,y,z) ~ f{t,x,y',z)\ ^ Kf^y \y ~ y'\ ; 
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2. for each {t,x,y,z,z') e [0, T] x R'' x R x M.'^'"^ x R^^^ 

\f{t, X, y, z) - fit, X, y, z')\ ^ (c + l{\zt + |z'|')) \z - z'\ ; 

3. for each {t, x, x' , y, z) € [0, T] x R'' x R'' x R x Ri^'^, 

\g{x) - g{x')\ ^ (C + |(|x|^/' + \x'\^^')) \x ~ x'\ ; 



4. 



a + Tl3 < 



Sometimes we will also consider stronger assumption. 

Assumption (B.2). Let / : [0, T] x R'' x M x R^^^ ^ R and 5 : M'^ ^ R be continuous functions and let us 
assume moreover that there exist five constants, Z^l, 0^r<j, a^O, /3^0, 7^0 and A'/.j, ^ such that: 

1. for each {t,x,y,y',z) e [0, T] x R^ x R x R x R^""^, 

\f{t, X, y, z) - f{t, X, y', z)\ ^ Kf^y \y-y'\; 

2. for each {t,x,y,z,z') G [0, T] x R"* x R x R^'^'' x Ri^"*, 

\fit,x,y,z)- f{t,x,y,z')\ ^ (c + ^(|;^|' + |z'|')) N - ^'1 ; 

3. for each {t, x, x' , y, z) G [0, T] x R'' x R'' x R x Ri^'^, 

X, y, z) ~ fit, x', y,z)\^(^C+ ^{\x\' + jxT)) \x - x'\ , 

\g{x)-gix')\^[c+^i\x\^ + \x'n) \x~x'\. 

Remark 2.1 Assumption (B.2) implies assumption (B.l). Moreover, the quadratic case corresponds to I = 1. 

Proposition 2.2 We assume that assumptions (F.l) and (B.2) hold. There exists a solution (Y, Z) of the Markovian 
BSDE (EU in S'^ x M"^ such that, 

\zt\^c{i + \xtn. 

Moreover, this solution is unique amongst solutions (Y, Z) such that 
• there exists ry > such that 



E 



< +00. 



Proof of the proposition First of all, let us remark that if we have a solution (Y, Z) such that 

\Zt\i^C{i + W), Vie[o,r], 

then, for all c > 0, 



E 



< +00 



(2.1) 



because 2lr < 2 (see e.g. part 5 in ||6l). Now, let us start by the uniqueness result. We consider two solutions 
iY\Z^) and iY^,Z^) such that g 52^ ^ c{l + and there exists 77 > such that 



E 



< +00. 
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We define Y := — Y^ and Z := — Z^. By considering the difference of the two BSDEs, the classical 
linearization method gives us 



YsUs + Z,Vsds~- [ ZsdWs, 
t Jt 



that is to say 

Yt = -J^ e^*'^-''"Z,{dWs-Vsds) 
where {U, V) takes value in M x M'' and 

By applying Young's inequahty, Holder's inequality and (12.1b . we have 



(2.2) 



Vs\^c + l{\zlf + \z',f). 



E 



^ E 

^ CE 

^ CE 

< +0O, 



:j;^(c+c\zl\"+{i+v)4\zT^'i- 



,Cl,^\zl\''ds^l^4j;f\zi\''ds 



1/p 



with q ^ (1/2 + r/)(l/2 + This estimate shows us that Novikov's condition is fulfilled and so we are 

able to use Girsanov's Theorem in (I2.2l i that gives us directly that Y = Q. Then it is standard to show that Z = Q. 
Finally we obtain the uniqueness result. 

Now, let us show the existence result. Firstly we will approximate our Markovian BSDE by another one. Let 
(yM the solution of the BSDE 

y,^ = sm^Xt) + /m(s, X„ r/^ )d,s - Z^dW,, (2.3) 

with gi\j ^ go p]\j and fuj = ./ ( • , Pm (•):•:■) where pM is a smooth modification of the projection on the centered 
euclidean ball of radius M such that \pm\ ^ M, \Vpm\ ^ 1 and pm{x) = x when |a;| ^ M — 1. It is now 
easy to see that gM and /m are Lipschitz functions with respect to x. Theorem 3.1 in 1251 gives us that Z^^ 
is bounded by a constant Aq that depends on M in the quadratic case. In fact this result stays true in our more 
general framework. More precisely we have this proposition that we will show in the appendix. 



pi xd 



and 



Proposition 2.3 We assume that (F.l) holds. We also assume that / : [0, T] X R'^ X R 

5 : R'' — > R are continuous functions such that: 

• g is K g-Lipschitz, 

• / is Kf x-Lipschitz with respect to x, K f y-Lipschitz with respect to y and locally Lipschitz with respect to 



z: there exists an increasing function ip 



such that for each {t, x, y, z, z') G [0, T] x 



\fit,x,y,z) - f{t,x,y,z')\ < C{1 + ^(|z|) + ^(|z'|)) \z - z'\ . 

Then, there exists a unique solution {Y, Z) to the BSDE ( 17.21 ) in 5^ x A^^ such that Z is bounded. Moreover, 
we have 

Thanks to this lemma we know that there exists a unique solution to the BSDE ( 12.3b (in the 

appropriate space) and is bounded by a constant Aq that depends on M. Moreover, /a/ is a Lipschitz function 
with respect to z and BSDE (12.3b is a classical Lipschitz BSDE. Now we will show the following lemma. 

Lemma 2.4 We have, 

|zf I ^a^ + bmW 

with {An, Bn)n&i defined by recursion: Bq — 0, Aq defined before, 

Bn+l = C, 

A„+i=C(1+A;'), 
where C is a constant that does not depend on M. 
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Proof of the lemma Let us prove the result by recursion. For n = we have akeady shown the resuh. Let us 
assume that the resuh is true for some n S N* and let us show that it stays true for n + 1. Firstly we assume that 
for all t e [0, T], b{t, .), g and /(i, ., ., .) are sufficiently differentiable. Then X and {Y^ , Z^') are differentiable 
with respect to x (see e.g. |l4l), we have 



+ £ v,/m(s, X,, y/^ )vx, + v,/m(5, x„ y/^ )vy/^ + v,/m(5, x„ y/^ zf')\/z^'ds, 

and Zl"^ = Vy/^(VXt)-V(t) a.s.. Since < Aq, we have 

|v,/M(s,x„y*^zf )| ^ c(i + |zf I') ^ c 

and so we are allowed to apply Girsanov's Theorem: Wt := Wt - V^/a/(s, X^, y/^ ^i^O^^s 

is a Brownian 

motion under a probability Q*^. We obtain 



vy 



and finally 



oo t 



alXrl' + pJ^ \X,\''ds 



because VXs{VXt) ^ is bounded by e-^'*^^ *^ Let us come back to the SDE: we have 

Xs^Xt+ I b{u,Xu)du+ f a{u)dWu+ f ct{u)V J{u, Xu,Y^' , Z^')du, 



\Xs\ ^ \Xt\+C+ / Kb\Xu\du 



a{u)dWu 



'^locT / |A„ + B„|X„r|'du, (2.5) 



E^ [\X,\] ^\Xt\+C + Kb J E^ du + CA'„ + CEY 

thanks to the recursion assumption. Young's inequality gives us 

Bl^\Xu[^ 

p q 



BllXul'^'du 



with l/p+l/q = l and rlq = 1 (let us recall that rl < 1). Thus, we obtain 



[\Xs\]^\Xt\ + C + Kb 
Gronwall's Lemma gives us 



]Xu\] du + CA'^ + CB!p + CEf 



\X„ \ du 



E? [\X, \]^C + CAi + CB^f + C\Xt 



and so 



E?" [\Xsn ^ (e?" [\xa]) cAii + cb:^p + c ix,r 

because r < 1. By introducing this inequality into ( 12.41 ) we obtain 

Izf U C + CAl^^ + CBll^P + C iX^r . 



Finally, we set 



Bn+l — C, 



(2.4) 
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and, thanks to the recursion assumption, we can take 

When h, g and / are not differentiable, we can prove the Lemma by a standard approximation and stabiUty results 
for Lipschitz BSDEs. □ 
Since Ir < 1, the recursion function that defines the sequence (A„)„^o is a contractor function and so A„ — > 
Aao when n — )- +00, with A^o that does not depend on M. Finally, 



\Z''\^A^+C\XX . 



(2.6) 



Now, we want to come back to the initial BSDE (11.2b . We will show that {Y", Z^)„eN is a Cauchy sequence 
in the space x Ai^. We have, thanks to the classical linearization method. 



that is to say 



with p,q e n, ^ and \Vp^'^\ ^ ^(1 + IZff + \Zp+'}\ ). Thanks to (1221), Novikov's condition is 

fulfilled and so we are able to apply Girsanov's Theorem: 



[f,+,is, X,,Yr\ ZP+^) Uis, X,,Yr\ ZP+^)] ds 



' Ef " [\u+,is,x,,YP+',zp+'^) - /p(s,x„r/+^zr?)|] 



(i + ixTr+')iix,i>p-i 



1 + \Xt\ 



2r+2 



't ['^\Xt\>p-1 



1/2 



(i + ix,r-+i)iix 



>p-i 



ds 



C / E> 



1 + 



E^ 



\Xs\>P~l\ 



,1/2 



ds 



< ceY 



1 + \Xt\ 



2r+2 



1/2 Ef " WXrW"^ 



C E> 



1 + 



|2r+2 



pl/2 

1/2 Ef' 



1/2 



By the same calculus than in the proof of Lemma (12.4b using the a priori estimate on ZP'"^ and ZP, we are able to 
show that 

Ef" [\X,n^C{l + \Xtr) 



for all a ^ 1 with a constant C that depends on a but does not depend on p and q. Finally 
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By applying Ito's formula to the process — Fp| and using the same calculus, it is rather standard to show 

that (Z")„gN is a Cauchy sequence in A4^. Finally, it is easy to check that (F", Z") "^^°°> (Y, Z) the solution 
of the initial BSDE (11.2b and our estimate ( 12.61 ) on stays true for Z: 



\Zt\ ^ Aoo + C \Xt\'' . 



□ 



Theorem 2.5 We assume that assumptions (F.l) and (B.l) hold. There exists a solution (Y, Z) of the Markovian 
BSDE such that 

s; C + e(^'(i+i/')+^/-)(^-*)(a + /3T)|a|^ei/'|Xt|i/', VtG [0,r]. 
Moreover, this solution is unique amongst solutions (y, Z) such that 

• there exists rj > such that 



E 



< +00. 



Proof of the theorem We will mimic the proof of Proposition (12.2b . Let us start by the uniqueness result. We 
consider two solutions Z^) and (y^, Z^) such that F\ € and 



E 



E 



^i2+v) 4 So \zi\"ds 



< +00. 



As in the proof of Proposition (12.2b we consider the BSDE satisfied by processes Y and Z and we introduce the 
two processes U and V. Now we just have to show that Novikov's condition stays fulfilled: by applying Young's 
inequality and Holder's inequality we have 



E 



^ E 
< CE 

^ CE 

< 



:j;^(C+{2+v)4\zT + i2+n)4\zT)ds 



1/2 



E 



1/2 



oo. 



The remaining of the uniqueness proof is unchanged. Now, let us show the existence result. We will consider 
again the solution {Y^^ , Z^^) of the BSDE ( 12.3b . We have already remark that Z^'^ is bounded by a constant 
Aq that depends on M. Now we want to obtain an estimate on Z^' that does not depend on AI by showing the 
following lemma. 



Lemma 2.6 We have, 



^A^{t) + B^{t)\Xt\^'\ 



with {An, Bn)n&i defined by recursion: Bq = 0, Aq defined before, 

Bn+i{t) - |a|^(a + /3T)e(^''(i+i/')+^/-)(^-*)e2'"'l-loc'^^^^W/', 
An+i{t) = kl^ {a + pT)e^K,(i+i/i)+K,„)(T-t)^2^-^w\^^TBl(t)/i (^(j ^ ^i-i/i |^|i/i ^i/'Ti/'A„(t)) 
where C is a constant that does not depend on M and t. 
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Proof of the lemma Let us prove the result by recursion. For n = we have akeady shown the resuh. Let us 
assume that the resuh is true for some n S N* and let us show that it stays true for n + 1. Firstly we assume that 
for all t e [0, T], b{t, .), g and f{t, ., ., .) are sufficiently differentiable. Then, by the same calculus than in the 
proof of Lemma|231 inequalities ( 12.4b and (12.51 ) stay true and we easily obtain 

Ef [\Xs\] ^ \Xt\ + C + K, j\T' [\X^\]du + 2''' |aL7r<(t) + 2'-i \a\^-fBl{mT' [ \Xu\du 
because t n- i?„(t) and t ^ An{i) are not increasing functions. Gronwall's Lemma gives us 

and so 

i/i 



^ i^f [\Xs\]) ^ (c + 2i-i/'|a|^S'/'Ti/'A„(0 + |Xt|^/')eK^^+2'-^kU7Bi(t))(T-t) 
because l/l < 1. By introducing this inequality into (I2.4l i we obtain 

Finally, we can take 
and, 

An+iit) = \<J\^ (a + /?r)e(^^(i+i/0+i^/,„)(T-t)g2'-i|<x|„-,TB;.(t)/i ^ ^i-i/i ^i/'Ti/'A„(t)) . 

When b, g and / are not differentiable, we can prove the lemma by a standard approximation and stability results 
for Lipschitz BSDEs. □ 
Now we want to study the behavior of the sequence (£?„ (i))„gN- Let us denote 

Ci(t) := \a\^ {a + pT)e^K,ii+m+K,,,)iT-t) ^i-i |^|^ ^j,^ 

Then, we have Bn+i{t) = Ci{t)e~^^'"^*\ It is easy to see that this sequence has a finite limit Boo{t) G if 
X > 0|a; = Ci{t)e~^ > is not empty and in this case we have 



Br 



,{t) = inf |a; > 0|a; = Ci(t)e^^' } 



Moreover, the set |a; ^ 0|a; = Ci{t)e | has only one, two or three elements. If |a; ^ 0\x = Ci{t)e ^ | 
has only one element denoted by x, necessarily we have 



f x = Ci{t)e^^' 
\ 1 = Ci(<)C2i'-^e't^* 

that gives us 

£ = Ci(i)ei/' 



Ciit) 



1 



Since we assume assumption (B.l).4, we have 



(eC2)i/'' 



so we conclude that the sequence has a finite limit Bao{t) < Ci(t)e^/' that does not depend on M. 

Now, let us see what happens to the sequence (A„(t))„gN. We can remark that 

A„+i(t) = Bn+i{t) (C + Cl"A^{t) 
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Since B„{t) Boo{t) < Ci(t)e^/' and Ci(i)e^/'C2^' < 1 due to assumption (B.l).4, a classical result for 
sequences defined by recursion gives us that (j4„(t))„gN has a finite limit Aoo(i) = — C'-Boc.(t)^^^ jj^^j. ^^^^ j^^j. 



depend on M because it does not depend on the initial value Ao{t). Finally, we have 



(2.7) 



To conclude we have to show that {Y", Z")„(=n is a Cauchy sequence in the space x A4^. As in the proof 
of Proposition l2.2l we consider the BSDE satisfied by processes yP+'J — Y'p and ^p+'J — and we introduce two 
processes t/^ ' and yP ''. Now we just have to show that Novikov's condition stays fulfilled: by applying Young's 
inequality and Holder's inequality we have 



E 



62 



^ E 
^ C] 



^ CE 

Thanks to estimate (I2.7l i we have, for all AI G N, 



1/2 



1/2 



(2.8) 
(2.9) 

(2.10) 



E 

< CE 
Since we have 



(2.11) 



sup ^ sup sup \X. 



a slight modification of the proof of Lemma 4. 1 in ifTSl gives us that the right term in inequality (12.111) is finite 
when 

which is true when 77 is small enough because we assume in assumption (B.l) that 

1 



a + T(3 < 



that is to say. 



1 1 

< 



2 V ^ I loo 2HLr22('-i) 2|a|Lr 

Finally Novikov's condition is fulfilled and we are allowed to use Girsanov's theorem. The remaining of the 
existence proof is unchanged. To conclude, we have to show that when we have a solution that verifies 

then, there exists 77 > such that 



E 



< +00. 



This claim is easy to prove because we have 



E 



^ CE 



cxp 



(l + ,,)re2UA'.(l + l/0+^/,„)T(^^^y)2;|^|2i g2y |;f^|2 



2 O^t^T 

and we have already shown that the right term in previous inequality is finite for 77 small enough. 



□ 
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Remark 2.7 In (B.l), our assumption 

a + TP< 



g^l/ig((l+/-i)/-fb+X/.JT |cr|^^/' TV 

is more restrictive than the one we can find in the article l\13]l for the quadratic case (i.e. 1 = 1), that is to says 

1 



a + r/3 < 



In this case, the assumption in l[13]l is more or less optimal because we need it to obtain a sufficient exponential 
moment for the terminal condition and the random part of the generator Let us also remark that in I113]I assump- 
tions are more general because they are about the growth of f and g instead of the growth of derivatives of f and 
9- 

Remark 2.8 With the same machinery it is possible to treat a little more general framework than the one of 
assumption (B.l): indeed it is possible to replace points 2 and 3 with 

2. for each {t,x,y,z,z') e [0, T] x R'^ x R x m.^""^ x Ri^'^, 

\f{t,x,y,z)-f{t,x,y,z')\ ^ (c + \x\ + ^Hzf + \z'f)) \z^z'\; 



3. for each (t, x, x' , y, z) G [0, T] 



X M'' X M'' X R X Ri^'^, 



\f{t,x,y,z)-f{t,x',y,z)\ ^ (^C + [3' \z\ + ^{\x\'/' + \x'\'/')^ \x~x'\, 

\g{x) - g{x')\ ^[C+ ^i\x\'/' + \xf/')) \x ~ x'\ ; 

and the point 4 with an Ad hoc assumption. We decided to do not deal with this little more general setting because 
the proof is already technical and we do not want to complicate it unnecessarily. 



3 Some results when a is random 

The main restriction in the previous part is about the function a that is assumed to be deterministic. In this section 
we will give some partial results when the SDE is given by 

Xt = x+ I h{s,Xs)ds+ f a{s,X,)dWs. (3.1) 
Jo Jo 

We will consider classical assumptions on this SDE. 

Assumption (F.2). Let 5 : [0, T] x R'' ^ R'* and ct : [0, T] x R'^ R'*^'^ be continuous functions and let us 
assume that there exist Kf, ^ 0, ^ and Ma- ^ such that: 

1. yt e [0,T], \b{t,0)\ C, 

2. Vt e [0, T], y{x, x') e R'* X R^ \b{t, x) - b{t, x')\ Kb\x - x'\, 

3. Vt e [0, T], Vx e R'^, \a{t, x)\ AU, 

4. Vt G [0, T], y{x, x') e R'* X R^ |cr(t, x) - a{t, x')\ ^Ka\x-x'\. 

Before giving our first result, let us point out why we are not able to use the same machinery than in our first 
part. When a is deterministic, the starting point is Proposition 12.31 where we show that Z is bounded under good 
assumptions. To prove this result we deeply use the fact that WX is bounded. Now, when a is not deterministic, 
WX is not necessarily bounded and so Proposition 2.3 does not necessarily remain. Finally, the first question to 
answer is; does the process Z remains bounded when g and / are Lipschitz with respect to x? 
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3.1 Boundedness of Z when T is small enough 

In this part we will give a partial answer to the previous question. 

Proposition 3.1 We assume that (F.2) holds. We also assume that f : [0, T] x R'^ x K x R^^'^ R and 
g : R'' — > R are continuous functions such that: 

• g is K g-Lipschitz, 

• / is Kf x-Lipschitz with respect to x, Kf y-Lipschitz with respect to y and locally Lipschitz with respect 
to z: there exists an increasing continuous function ip : R'^ — )■ such that for each (t^ x^y, z, z') € 
[0,r] X R'' X R X Ri^^ X Ri><^ 

|/(t, X, y, z) - fit, X, y, z')\ (/^/,, + + ip{\z'\)) \z - z'| . 

Then, for T small enough, there exists a unique solution (Y, Z) to the BSDE ( 17.21 ) in x such that Z 
is bounded. 

Proof of the proposition Once again, we will use a classical truncation argument (see e.g. the proof of Theorem 
4. 1 in |U2| ). Our truncation function pM is a smooth modification of the projection on the centered euclidean ball 
of radius M such that \pm\ < M, |VpA/| < 1 and pA/(a;) ^ x when |a;| < A/ - 1. We denote (F^, Z^^) the 
solution of the BSDE 



Y,^' ^ giXr) + fM{s,Xs,Yf,Zf')ds-J^ Z^ 



where /a/ := /(., ., ., pA/(-))- Now, this BSDE is also Lipschitz with respect to z. Firstly we assume that for all 
t G [0, T], b{t, .), g and f{t, ., ., .) are differentiable. Then X and {Y, Z) are differentiable with respect to x, we 
have 

Vr/^ = Vg{XT)VXT ~ / VZ^'dWs 

+ v,/a./(s, Xs, y/^ )vx, + v„/a/(s, x„ y/^ )vy/^ + v,/a/(s, x„ y/^ )vz 

and Z^'' = Vy/'^(VXt)"V(t,Xt) a.s.. Since V^/a/ is bounded by X/,^ + 2ip{M), we are allowed to apply 
Girsanov's Theorem: Wt := Wt — /g Vz/a/(s, y/'', Z^^)ds is a Brownian motion under the probability Q^^. 
We obtain 



and finally 



\Z,''\ < e^/-^A/, ( KgEf [|C/t|] Ef [\Us\]ds 



with {Us)t^s^T the solution to the SDE 

Us = Id+ f Vb{Xu)Uudu+ [ Y.V<7\Xu)UudW'^ 

= Id+ f Vb{X^)U^du+ f Y,^'^'iX-)UuidW^ + iVjMy{u,Xu,Y^',Z^')dv 

■li Ji i=l 



where the superscript ' denotes the i-th column. A classical estimate on [\Us\] gives us 
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For T small enough, the function 

has, at least, one positive fixed point. Let us denotes M the lowest positive fixed point of this function. Then 
Z^^ ^ M — 1 and {Y^^ , Z^^) is a solution to the initial BSDE. Uniqueness follows from the uniqueness result 
for Lipschitz BSDEs. □ 

Remark 3.2 

• In general, it is not possible to stick local solutions to obtain a solution (Y, Z) with Z bounded for all T. 

• The biggest T that allows the existence of a fixed point for the function 

strongly depends on Kg. So, it is not possible to treat the case of g and f locally Lipschitz with respect to x 
by using the same machinery than in the previous part. 



3.2 A simple example 

In this part we will see that when we consider a simple quadratic BSDE with an explicit solution, the process Z 
remains bounded when g and / are Lipschitz with respect to x. More precisely, we will consider the following 
quadratic BSDE; 

Yt = g{XT)+ [ ^-^ds- [ Z.dWs. (3.2) 



Proposition 3.3 Let us assume that (F.2) holds and that g is a K g-Lipschitz function. Then there exists a unique 
solution (y, Z) to the BSDE ( 13.21 ) in x Ai^ such that the process Z is bounded. 

Proof of tlie proposition It is well known that (I3.2l i can be explicitly solved with an exponential transform, also 
called Cole-Hopf transform in PDE theory. More precisely, we have 



Yt^logEt 



, Zt^e-'^'Zt 



where Z is given by the martingale representation theorem applied to the martingale (Ej [e^^^^^-*] )q^j^j,. Y is 
well defined because g is Lipschitz and for all C > 



< +0O, 



since <t is bounded (see e.g. part 5 in |6j|). The uniqueness is standard. As in previous proofs, we assume in a first 
time that g, b and <t are differentiable with respect to x. Then we have 



\Zt\ 



Et [Vg{XT)VXTiVXt)-^ 



-lpS(^T)l 



s; c 



Et [esi^T^ 
Et [|VXT(VXt)-i| e9(-^^)] 



-a{t,Xt) 



CEt 



Et [eaiXT)^ 
\VXTiVXt)-^ 



L/2 Et [e^siXT)^ 



1/2 



^ C 



Et [e^aiXT)^ 



1/2 



Et [e9(^r)] 

because Vg, a are bounded and {V XsC^ Xt)^^)t^s^T solve the SDE 

d 



Us^Id+ f Vb{u,Xu)Uudu + J2 I '^a\u,Xu)UudWl 

Jt ; 1 Jt 
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so Et 



random initial condition) 



is bounded. Let us denote X the solution of the ordinary differential equation (with a 



Xs = Xt + J b{u, Xu)du, ts^si^T. 



Since g is Lipschitz and Xt is J^t -measurable, we obtain 



\Zt\€C- 



2,C\ X'jp — X'j' 



1/2 



Let us estimate 



2C\Xt-Xt\ 



. We have 



-C —X.'j 



C- 



'2C \ X'j' — X'j' 



1/2 



- C Xt — X'J 



(3.3) 



\Xs- Xs\=Q + Kb I \Xu-Xu\du+ sup 



a{u,Xu)dWu 



and we deduce from Gronwall's lemma the inequality 

a{u,Xy,)dWu 

It follows from the Dambis-Dubins-Schwarz representation theorem that, for A ^ 0, 



sup Xs — Xs ^ C sup 

t<s<T t<s<T 



E 



exp ( A sup 

t<s<T 



a{u,Xu)dWu 



) 


J^t,Xt = xo 




cxp A sup 


J a{u,X'^^°)dWu 













s$ E 



sup e 

O^ssSlkllLT 



X\Ws 



< +00 



where {Xl'^°)t^s^T stands for the solution to the SDE ( 13. Il l that starts from xq at time t. We remark that the right 
term in the last inequality does not depend on xq so 



E. 



exp ( C sup 

t<s<T 



a{u, Xu)dWu 



is upper bounded. By the same type of argument we have that 

exp — C sup / a{u,Xu)dWu 
is lower bounded by a strictly positive constant. Finally ( 13.31 ) becomes 



\Zt\ ^ c- 



E. 



-C sup,^^^^\l^ a{u,X^)dW^\ 



and so Z is bounded. Finally, when g, b and cr are not differentiable, we can prove the result by a standard 
approximation. □ 

Remark 3.4 Thanks to this estimate on Z, it is possible to use the same machinery than in the previous section to 
show estimates on Z when g and f are locally Lipschitz with respect to x. This simple example is a good argument 
to postulate that Theorem \2.5\ or Proposition \2.2\ could stay true when we replace (F.l) by (F.2), at least in the 
quadratic case. 

3.3 The case of bounded terminal conditions 

In this part we will restrict our study to the quadratic case and we will assume that the terminal condition and the 
generator are bounded with respect to x. In this case we are able to obtain estimates on Z thanks to the additional 
tool of Bounded Mean Oscillation martingales (BMO martingales for short). We refer the reader to Il20l for the 
theory of BMO martingales and we just recall the properties that we will use in the sequel. Let = Jq (j^sdWs, 
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for t G [0, r], be a real square integrable martingale with respect to the Brownian filtration. Then $ is a BMO 
martingale if 

-,1/2 



I* 



\bMO^ sup E[($)t- Wrl^r]'/' = SUp E 

tg[o,t] Te[o,T] 



ids 



< +00, 



where the supremum is taken over all stopping times in [0, T] and ($) denotes the quadratic variation of $. In our 
case, the very important feature of BMO martingales is the following lemma: 

Lemma 3.5 Let ^ be a BMO martingale. Then we have: 

L The stochastic exponential 



£{^)t = ft = exp (^J^ (j>sdWs - i ^ 10.1' ds 



s$ t T, 



is a uniformly integrable martingale. 

2. Thanks to the reverse Holder inequality, there exists p > 1 such that £t G L'p. The maximal p with this 
property can be expressed in terms of the BMO norm of^. 

We will work under following assumptions on coefficients of SDE ( 15.61 1 and BSDE ( 11.21 ). 



Assumption (F.3). Let b : [0, T] 



and cr : [0, T] x 



pdxd 



be continuous functions and let us 



assume that there exist Ki, ^ 0, ^ 0, ^ and k G [0, 1] such that: 

1. vt € [o,r], |6(t,o)| c, 

2. Vt G [0, T], W{x, x') e M'* X |6(t, x) - b{t, x')\ s$ Kb \x - x'\, 

3. Vt € [0,T],yx e M'*, \cr{t,x)\ A/^(l + 

4. Vt G [0,T],y{x,x') e M'* X \a{t,x) ~ a{t,x')\ ^ K„ \x-x'\. 

Assumption (B.3), Let / : [0, T] x M'^ x M x Ri^"^ ^ M and 5 : M'* ^ M be continuous functions and let us 
assume moreover that there exist seven constants, r G IR+, a ^ 0, /? ^ 0, 7 ^ 0, K f^y ^ 0, Mf ^ and Alg ^ 
such that: 

1. for each {t,x,y,y',z) e [0, T] x M"^ x M x M x R^""^, 



\f{t, X, y, z) - f{t, X, y\ z)\ Kf^y \y-y'\\ 



2. for each {t,x,y,z,z') £ [0,T] x x R x M.^""^ x IRi^'*, 



\fit,x,y,z)~fit,x,y,z')\ ^ (c + ^(|2| + |z'|)) \z ~ z'\ 



3. for each (t, x, x', y, z) S [0, T] x M'' x E'^ x E x R^'"', 



|/(t, X, y, z) - fit, x', y,z)\^(^C+ ^{\x\' + jxT)) \x - x'\ , 

\g{x)-g{x')\^[c+^i\x\^ + \x'n) \x-x'\; 

4. for each (t, x, y, z) € [0, T] x E'' x E x R^""^, 

\f{t,x,y,z)\^Mf{l + \y\ + \z\'), 
\g{x)\ ^ My. 

Theorem 3.6 We assume that assumptions (F.3} and (B.3) hold. There exists a solution (Y, Z) of the Markovian 
BSDE in S'^ x and this solution is unique amongst solutions (Y, Z) E x such that Y is bounded. 
Moreover we have 

|ZtKC(i + |Xtr+'=), vtG[o,r], 

and 



BMO 



where the last constant C depends only on Mg, Mf and Kf^,, 
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Proof of the theorem For the existence and uniqueness result we refer the reader to 11211 123]| . The estimate for 
the BMO norm of Z is shown in H [l]. It just remains to prove the estimate on Z. As in previous proofs, we 
firstly assume that /, g, b and a are differentiable with respect to x. Then, according to |I4][T1, X and {Y, Z) are 
differentiable with respect to x, we have 

+ \/^f{s,Xs,Ys,Z,)\/Xs + \7yf{s,Xs,Ys,Zs)\7Ys + VJ{s,X„Y„Z,)\/Z,ds, 
and Zt = VYt{VXt)-^a{t,Xt) a.s.. Since Z^dWs is BMO and 

|V,/(s,X„n,Z,)| sc: C{i + \Zs\) 

then S/zf{s, Xg, Yg, Zs)dWs is BMO and we are allowed to apply Girsanov's Theorem thanks to Lemma 
Wt := Wt — /g ^zfis, Xg, Yg, Zg)ds is a Brownian motion under the probability 

Q^£(^J^Wzf{s,Xg,Yg,Zg)dWg^ P. 



We obtain 



ei: v„/(«,x„,Y„,z„)d«y^^(^^ j^^^Yg, Zg)VXgds 



1/2 



1/2 



sup |VX,(VXt)-M 

t<s<T 



-,1/2 



(3.4) 



and then it comes 

l^tKC^l+E? [iXtI''' ■ + / \Xgr ' rfslE^ 

by using assumptions (F.3), (B.3) and Cauchy-Schwarz's inequality. Let us denote 

ft,T:=exp|^^ VJ{s,Xg,Yg,Zg)dWg~^J^ \V Jis, Xg,Yg, Zg)\^ ds^ . 

Thanks to Lemma [3751 there exists p > 1 (that does not depend on t) such that Et[£^j,] < +oo. But, by using 
Holder's inequality and classical estimates on SDEs we have 



E^ 



\Xg 



i2rg 



and 



E> 



sup \wXg{wxty 

t<s<T 



sup \vxg{wxty 

t<s<T 



-1 |2g 



1/9 



By putting th two last inequalities into (13.41 ) we obtain the result. Finally, when b, g and / are not differentiable, 
we can prove the result by standard approximations and stability results for quadratic BSDEs (see e.g. ifTSl ). □ 

4 Application to quadratic and superquadratic PDEs 

In this section we give an application of previous results concerning BSDEs to semilinear PDEs which have a 
quadratic or superquadratic growth with respect to the gradient of the solution. We will restrict our study to 
deterministic functions cr. Let us consider the following semilinear PDE 



dtu{t, x) + £u{t, x) + fit, X, u{t, a;),* Vu{t, x)a{t)) =0, x eM.'^,te [0, T], 
u{T, .) = g, 



(4.1) 
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where C is the infinitesimal generator of the diffusion X*'^ solution to the SDE 



The nonlinear Feynman-Kac formula consists in proving that the function defined by the formula 

V(t, x) e [0, T] X R'', u{t, x) y/'"^ 



(4.2) 



where, for each {to,xo) e [0,r] x R'^, (Y^o.xo ^ ^to,xo-^j stands for the solution given by Theorem 1231 to the 
following BSDE 



is a solution, at least a viscosity solution, to the PDE (14.1b . Firstly, let us study the growth and the continuity of 
this function. 



Proposition 4.1 Let assumptions (F.l) and (B.l) hold. The function u defined by i4.2\) has a polynomial growth 
and is a continuous function. More precisely we have, y{t, t' , x, x') G [0, T]^ x x W^, 

\u{Ux)\ ^ + 

\u{t,x) - u{t', x')\ ^ C(l + \x\'^' + \x - x'\ + C(l + \x\'+'^' + |xf \t - t'\^^^ . 

Proof of the proposition To show the first point, it is sufficient to prove the estimate 



E 


sup 











2(1 + 1//) 



By a very classical method we can easily show the estimate 

e[ sup sc CE \g{X'^^)\^+ r \f{s,Xl^^,0,Zln 

t^s^T J Jt 

Since ^ C(l + we obtain, by using the growth of g and / and classical estimates on SDEs, 

1+ sup Ix^l'^'^'/'^^ 

t<s<T 



(4.3) 



E 


sup 


sc CE 









Now, let us show the second part of the proposition. By a symmetry argument we are allowed to suppose that 
t' > t. Then 

u{t, x) ~ u{t', x') = E [y/'^ - r,*'^] + E 



Cauchy-Schwarz's inequality and growth assumptions on / and g give us 



E^* f{s,Xl'\Y:'\Zlnds 



|E [y/'" - y,*'"] f = 

s$ \t-t'\E 
s$ C\t-t'\E 

Thanks to a priori estimate on Z, a classical estimate on SDEs and (14.31 ), we obtain 

|E [r/'" - y,*'"] f ^ c |t - t'l (1 + 



\fis,Xl'^X-^,Zl'^fds 

1+ sup (|x,|'(i+i/') + |y/ + |^.,|" 

f<s<T ^ 
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Now we will study the term 1 



We have, thanks to the classical linearization method, 



with 



^xx 



^ A'/_j^ and 



yx.x 



7t'x' 



). Since Novikov's condition is fulfilled, we are able 



to apply Girsanov's Theorem. We obtain, by using the fact that / and g are locally Lipschitz, 



Iff -f+z 



/(s, x*-^ y,*'^ z*^-) - fis, x^'^'X'^, z*f') 



ds 



sup ( 1 + \X 

t'<s<T 



t.x\Ul 



X 



t'x' 



l/l 



\rt.X \rt .X 



C 1 + 



sup \X 

t'<s<T 



■11/2 






2/r 


l/2\ 


+ Ef' ^ 


sup 



















\rt'X vt' .x' 


2" 


sup 











Let us recall that 
to obtain finally 



yx.x 



^ c 



1/2 



Once again we are able to use classical methods on SDEs 



\^t,X ~trt ,X 



Classical estimates on SDEs allow us to conclude: 



E 



\^t,X xz-t ^x 



□ 



Proposition 4.2 Lef assumptions (F.l) and (B.l) hold. The function u defined by A4.2i is a viscosity solution to 
the PDE KB- 



Since we are able to use Girsanov's transformation in the BSDE, we have a comparison result. Moreover, Propo- 
sition l4.ll gives us that u is a continuous function. So the proof of the proposition is totally standard: for example, 
it can be easily adapted from the proof of Theorem 4.2 in lfT4l . 



5 Time approximation of quadratic and superquadratic Markovian BS- 
DEs 

5.1 approximation of the initial BSDE by a Lipschitz one 

In a first time, we will consider the deterministic case for the function a and we will approximate the solution 
(y, Z) of the BSDE (O by (F^^, Z*^) the solution of the BSDE (jSj]). The aim of the following proposition is 
to study the approximation error given by: 











ei(Af) :=E 


sup iFt-r/^f 


+ E 










Jo 



Proposition 5.1 If we assume that assumptions (F.l) and (B.l) hold, then there exists A > such that 

ei(M) < Ce-^^^\ 
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Proof of the proposition Let us define processes SY :=¥ — and SZ := Z — Z'^' . We have 



5Yt = giXr) - 9m{Xt) + / /(s, X„ F,, Z,) - X„ F/^ , )ds 



The classical linearization method gives us 



5Yt ^6g+ / J/,, + SYsU^' + SZ.Vf' ds 



SZsdWs, 



(5.2) 



with 

5g:=giXT)-9M{XT), 5f, := f{s, X,,,Ys, Z,) ~ fuis, X,,Y,, Z,,), 
(jjM ^ yM-^ ^jjjj y^j^g in M X and 



We can easily show that Novikov's condition is fulfilled for V'^^ by doing the same calculus than for V^^'' in 
the proof of Theorem l2.5l (inequalities ( 12.81 ) to ( 12.10b ). So, we are allowed to apply Girsanov's theorem: Wt := 
Wt — /q V^^ ds is a Brownian motion under the probability Q^^. Thus, by applying Cauchy-Schwarz's inequality 
and Markov's inequality we obtain 



5Yt 



\SYt\ ^ cw!^ 



(1 + |Xt|'+'/')1|xh^a/ 



C 1 + E' 



\Xn 



|2(l + l/i) 



1/2 it^r 



{l + \Xs\'+'^%x^\^Mds 

ol 1/2 



+C / 1 + 



IX 



,2(1 + 1/0 



1/2 !i^t 



1/2 



-ds. 



Then we use the following lemma that we will prove in the appendix. 
Lemma 5.2 We assume that assumptions (F.I ) and (B.l) hold. We have 
• Va e [l,+oo[, 3C > 0, 



sup iXsl" 

.t<s€T 



^C{i + \Xtn, Vte[0,T], 



(5.3) 



3C>0, 3fi>0, V^iG [0,/i[, 



sup e 

t<s<T 



Now i5.3\ becomes. 



By using Cauchy-Schwarz's inequality, we obtain for all p ^ 1 and for all < A < A with A small enough, 

C 



E 



sup \SYtf 



c 

ppXAP 



-E 



(1+ sup |Xt|^(^+^/'))e'^f^™P»^*<^'^* 



1 + E 



sup \Xt 



2p(l + l/0 



l/2^ 



E 



=CpAsupo^j^y|Xt| 



1/2 
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Let us remark that C depends on A but does not depend on A. By using classical results about SDEs (see e.g. the 
beginning of part 5 in jgl) we have, for all p ^ 1, 



E 



sup \Xt 



2p(l + l/i) 



1/2 



< +00, 



and, for A small enough. 
Finally we obtain that 



1/2 



E 



sup \6Ytf 

OsSt<T 



< +00. 



C 



(5.4) 



To study the error on Z we come back to i5.2i and we apply Ito's formula: 



\5Yt\^+( \SZ,fds = \5YT\''+ [ 26Ys{Sfs + SYsU^' + 5Z,Vl')ds- [ 2SYsSZ,dWs 
Jt Jt Jt 



We obtain by applying Cauchy-Schwarz's inequality 



E 



\5ZtYdt 





















'\5Yt\^' 


+ 2E 


/ SYtSftdt 


+ 2E 


r \5Yt\'U^'dt 


+ 2E 


[i 






Jo 




Jo 







s$ E 



< (l + 2i^/,,)E 



SYtSZtV.'^'dt 



sup \SYt\ 

0<t<T 



2yi/2E 



-,1/2 



sup \SYt\ 



E 



1/2 







1/2 




1/4 


+2E 


/ \SZtfdt 


E 


sup \SYtf 


E 




Jo 









sup 



sup \Sft\ 

1/4 



.0<t<T 



Thanks to inequalities 
and 

it is easy to see that 









^C{1 


+ \x, 


1/') 


\SVt^'\ ^ c 






)^C(1 


E 


sup 


\S.ftf 


+ E 


sup 

Os^t^T 





s; c 



with C that does not depend on AI. Then, by applying iSAi and the inequality 2ab ^ ^ + 26^ we have 

-I 1/2 



E 



ISZtl'dt 



c 



c 



2E 



\SZt\'dt 



c 



Finally we obtain 



C 1^ 
^ ;^+2^ 



\SZt\Ut 



c 



(5.5) 

Jo J e'""" 

To conclude, (15.4b and ( I5.5l l give us the result. □ 
Now we want to obtain the same type of estimate in the quadratic case when a is random. Since cr is not 
necessarily bounded, Z could be unbounded even if g and / are Lipschitz functions with respect to x. So, we will 
approximate the solution (F, Z) of the BSDE (O by (F^^, Z*^) the solution of the BSDE 



= (Xt)) + I f{s,p^j„.^^,-.{Xs),Yf\PM{Z':'))d. 



M 



~zM\ 



Zl'dW, 



(5.6) 
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where pM is a smooth modification of the projection on the centered euclidean ball of radius M such that \pM \ ^ 
M, |V/9j\/| ^ 1 and pm{x) = x when |a;| ^ M — 1. The aim of the following proposition is to study the 
approximation error given by: 



ei(Af) :=E 


sup 


+ E 


r \zt-zt'\^dt 


(5.7) 








Jo 





Proposition 5.3 If we assume that assumptions (F.3), (B.3) hold and 2k ^ 1 — r, then there exists A > such 
that 

ei(Af) ^ Ce~^'''\ 

If moreover 2k < 1 — r, then there exist A > and £ > such that 

ei(M) <Ce-^*^'^\ 

Proof of the proposition Thanks to BMO tool, we have a comparison theorem that gives us an estimate for ei . 
Indeed, we can apply Lemma 3.2 in ifTSl : there exists q> 1 such that 



+CE 



|2q- 



1/9 



1/9 



(5.8) 



Assumptions (F.3), (B.3) and the estimate on Z give us 

+c(i + ix,r+'')i|^,.+». 



,r+''^M/c-i 



< c(i+ sup |x,r+'')i 

0^s<T 



and 



ir-t-K. 



< c(i+ sup )isup„<,<^|x,r+'<>j\//c-i- 



By using Holder's inequality and the fact that, for all p > 0, 



E 



sup \X,f 

0<s<T 



< +00, 



becomes 



ei{M) CP sup |X,r+'' ^ M/C - 1 

with q' > 1. To conclude, we will use the following lemma that will be proved in the appendix. 
Lenuna 5.4 We assume that (F3) holds. There exists A > such that 

i2(1-k) 



(5.9) 



E 



exp A sup \Xt\ 



< +00. 



Since we have assume that r + k ^ 1 — k, Markov's inequality and previous lemma give us, for M big enough, 



sup |x,r+" ^ M/c - 1 ) s=: rvv^^^^.. 

0<s<T 



=A(A//C-1)2 
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Then, the first part of the Lemma is obtained by putting this inequaUty in the estimate ( 15.91 1. When r + k < 1 ~ k, 
we denote e := — 1 > and we obtain, for M big enough. 



sup \Xt 



I r+K 



> M/C - 1 s; 



pA(M/C-l)2(i+^) 



c 



„AM2(i+6) 



Finally, the last part of the Lemma is proved by putting this inequality in the estimate i5.9i . 
Remark 5.5 When 2k > 1 — r it is possible to show with the same proof that there exists A > such that 

C 



□ 



ei 



exp 



When K = 1, it is also possible to recover the result obtained by Imkeller and dos Reis in l[18]l ( they assume in 
addition that r = 0): for all fc e N, there exists C > such that 

Let us remark that our result is more precise than the one of l\I8]l and our proof is more simple since we do not 
have to study the second order Malliavin differentiability of the BSDE. 



5.2 time approximation of the BSDE 

In a second time, we will approximate our modified BSDE by a discrete time one. We denote the time step by 
h ~ — and [tk ~ kh)a^k^n stand for the discretization times. One needs to approximate X by a Markov chain 
X" which can be simulated. For example, we will consider the classical Euler scheme given by 



XT 







X^^ + hb{tk,Xl) + a{tk,Xl){Wt,^, -Wt,), sc: fc n - L 



We denote (F^-", Z*^'") (resp. (y^^^", Z*^-")) our time approximation of (y^^ Z^O (resp. Z^^))- These 

couples are obtained by the classical explicit dynamic programming equation: 



M.i 



Ma 



s$ A: s$ n, 



^n, 



and 



Y 



M.n 
7M,n 



= .g(PM(.-+»)-(^t"))' 



fc n, 

yM,n 



In a classical framework, there is already results about the speed of convergence of BSDE time approximation. 
Let us precise the classical result shown by Il3l l29ll22l . 

Proposition 5.6 Let us assume that assumption (F.l) or (F.3) holds. We also assume that 

• g is K g-Lipschitz, 

• / is K f x-Lipschitz with respect to x, Lipschitz with respect to y and Kf _z-Lipschitz with respect to z. 

We denote (K", Z") the time discretization of {Y, Z) given by the classical explicit dynamic programming equa- 
tion. The error discretization is given by 



e(n) 



sup E 



Y" ~Y, 



n-l 

-El 

fc=0 



Zf dt 



Then, there exists a constant C that does not depend on Kg, Kf^ and Kfz such that 



e{n) s; Ce'^'^f 



l + K 



'2 



-E 











/ \Zt\'dt 


+ E 


sup |yt|' 




Jo 
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This proposition will be proved in the appendix. Now, the aim of this section is to study errors of discretization 
e(M, n) and e{M, n) given by 



EE 



tk + l 



7M.I 



dt 



and 



e(M, n) :— sup E 



e{M, n) := sup I 



Theorem 5.7 We assume that assumption (F.I ) holds and {Y , Z ) is the solution ofBSDE ( 12.31 
• Let assumption (B.2) holds. We have 

C CeC*^"' 







2' 


n-1 






2 


















fc=0 


Jtk 







/n particular, for all 1 < p < (rl)^^, if we take M = (logn)^/^ then e(M, n) — o{h^^^) for all e > 0. 
• Let assumption (B.l) holds. We have 



z{M,n) < 



C C, 



In particular, if we take M = -Qp^^/^ogn then e(M, n) = o{h'^i+^). 

Proof of the theorem It is easy to see that 

e{M,n) 2(ei(Af) + e2(A/,n)) 

with ei(A/) defined by (ISTTT i and 



e2(M, n) := sup E 



it. f 



The error ei{M) is already studied in Proposition l5.ll Concerning the error 62 (Af, n) let us remark that (F^^, Z'^^) 
is the solution of a BSDE with Lipschitz coefficients: indeed, and are Lipschitz with respect to x and y, 
Jm is locally Lipschitz with respect to z and Proposition 12 . 3 1 gives us that Z^^ is bounded. Thus, we are allowed 
to apply Proposition l5.6l 



•^0 



S^\Z^^\'dt] +E[supo^,^T|y/'|']) 



e2(M,n) 

with A'gjj^ the Lipschitz constant of gM, and Kf^,^x, the Lipschitz constants of fuj with respect to x and 

z. Estimations on given by Proposition 12.21 and Theorem 12 . 5 1 show us that E \ dt 



is bounded by 



a bound that does not depend on AI. Thanks to Ito's formula applied to e^^M"* li^i^^T and estimations on Z 



zM 



given by Proposition 12.21 and Theorem 12. 5 1 it is also possible to show that ] 
by a bound that does not depend on M. Thus, we have 



A/ I 



is also bounded 



e2(M,n) ^ 



Under assumptions (F.l) and (B.2) we have, thanks to Proposition l2.3l Kg^^, ^ C(l + il/'), Kf^^^x ^ C(l + Af) 
and s$ C(l + (C(l + Af))') s$ C(l + A/'^')- Finally, we obtain 



e2{M,n) ^ 
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Under assumptions (F.l) and (B.l) we have, thanks to Proposition 1231 Kg,^, C(l + M^/'), i^/^^.a; ^ C(l + 
Afi/') and Kf^,^^ ^ C(l + (C(l + Mi/'))0 «^ C'll + ^^)- Finally, we obtain 



e2(M,n) 



□ 



Remark 5.8 Since a is a deterministic function, Euler and Milshtein schemes are equal, so the discretization 
error on X is better In this situation, authors of ^I5]l show that the discretization error for the BSDE is on the 
same order than the discretization error for the SDE if we assume extra smoothness assumptions on b, a, g and f. 
More precisely, we could obtain the better estimate 



e2{M,n) ^ —. 

Theorem 5.9 We assume that assumptions ( F.3 ) and (B.3) hold and , Z^^ ) is the solution of BSDE 
• If we assume that 2k < 1 — r then there exists rj such that 

C Ce^^'" 



e(M, n) < 



In particular, for all (2 + 77) ^ < p < 1/2, if we take M = (logri)^ then e{M, n) = o{h^ ^) for all £ > 0. 
• If we assume that 2k = 1 — r then we have 

, C Ce^-^^" 
e{M,n) -^-^ + _— . 

In particular, if we take M = V^ogn then e{M, n) = o(h'^i+'^2 ). 

Proof of the theorem It is easy to see that 

e{M, n) ^ 2{ei{M) + 62 (M, n)) 

with ei(M) defined by (l5Jl l and 



62 (M, n) := sup E 

0<k<n 



Z. 



M 



dt 



The error ei{M) is already studied in Proposition l5.3l Concerning the error 62 (Af, n) let us remark that (?^^, Z*^) 
is the solution of a BSDE with Lipschitz coefficients: indeed, (.)) and /(., (.), ., pm{-)) 

are Lipschitz with respect to x, y and z. Thus, we are allowed to apply Proposition l5.6l 



K 



E 



/o \zn dt 



e2(M,n) s; 

li 

with Kg the Lipschitz constant of ,g(p^y^(i_^,-i (.)), and if Kf^^ the Lipschitz constants of /(., pj^^i^^)-! (.),., pm{ 



with respect to x and z. Classical estimates on solutions of quadratic BSDEs show us that E 



Ml 



andE 



Ml 



dt 



are bounded by a bound that does not depend on M. Thus, we have 



Under assumption (B.3) we have. Kg ^ C(l + A/''^!"'')), A"/,^ C(l + Af''(i-'=)) and if/,^ C(l + A/). 
Finally, we obtain 

(7gCM= 

e2(Af,n) s; . 

n 

□ 

Remark 5.10 When 2k > 1 — r, f/ie error estimate for ei(Af) given in remark \53\ is not sufficient to obtain a 
"good" speed of convergence: the estimate on e(Af, n) becomes, for M well chosen, 

eiM,n) < ^ 

(log nj'^ 

/or fl// k G N*. r/i/i phenomenon is already explained in introductions of articles 
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A Appendix 



A . 1 Proof of Proposition 12.31 

To show the resuh we will use a classical truncation argument (see e.g. the proof of Theorem 4.1 in [12]). Our 
truncation function is the projection on the centered euclidean ball of radius N in E^^''. We denote (F^, Z^) 
the solution of the BSDE 

Y/" = giXr) + r fis, X^^Y!" , pM[Z^))ds - C Zf dWs. 



Now, this BSDE is also Lipschitz with respect to z. By the same calculus than in the proof of Theorem 3.1 in ll26l 
we can show that Z^ is bounded by 

\Z^\^e'^'^^+^'-^^^\a\^{K, + TKf,^). 

This bound does not depend on N so pn{Z^) = for N big enough. Then a uniqueness result for BSDEs 
with Lipschitz coefficients gives us that (Y, Z) = (y", Z^) and the result is proved. □ 



A.2 Proof of Lemma 15.21 

Thanks to the estimate on Z of Theorem l2.5l we easily show 



Xs = + J h{u,Xu)du + J a{u)[dWu + Vfdu] 



\X,\ ^ \Xt\+C + C I \XJ du- 



ff {u)dWu 



sup \Xr\ ^ \Xt\+C + C / sup \Xr\du+ sup 



a{u)dWu 



and we deduce from Gronwall's lemma the inequality 

(7{u)dWu 

The first part of the lemma is easily proved thanks to the previous inequality. Moreover, we also have 



sup \Xr\ ^ C 1 + \Xt\ + sup 

t<r<s V t<r<T 









sup exp Cp 


J aiu)dWu 


■)] 









It follows from the Dambis-Dubins-Schwarz representation theorem and Doob's maximal inequality that 



sup exp Cp 

t<s<T \ 



which is a finite constant ifC/.j|(T|^T<l/2. 











J a{u)dWu 


11 


s$ E 


sup 











^ 41 



□ 



A.3 Proof of Lemma 15.41 

Let us consider the process 



Yf, :-- 



Ito's formula gives us 



Y 



Yo + {l- fi) 



f ds + \ f trace {V^F{X,)a{s, X,Ya{s, X,)) ds 

Jo (1 + \XJ ) 2 Z Jq 



lo {i + ix^n 

+ (1 - f -^^^^^^dw. 

Yo+ [ bis,Xs)ds+ I d{s,X,)dWs. 
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with IctI ^ C and 



b{t,x) C(l + Then, we have 



sup \Yt\ |yol + 

0<t<u 



sup 

O^t^s 



ds + sup 



< \Yo\+CT + C [ sup |Xf|^ "^5+ sup / a{s,Xs)dW, 

< |yo|+CT + c/ sup \Yt\ds+ sup / CT(s,X,)dM^, 

Jo O^t^s O^t^T Jo 



and we deduce from Gronwall's lemma the inequality 



sup \Yt\ s^C [1+ sup 



O^i^T 



0^t<T 



f a{s,Xs)dWs 

Jo 



Since a is bounded, we are able to fit the end of Lemma lSSj to show that there exists A > such that 

|2 



E 



exp ( A sup \Yt\ 

0<t<T 



< +00. 



Since \Xt\^^ '^^ ^ |Yi|, the proof is complete. □ 

A.4 Proof of Proposition 15.61 

It is already proved in ||3]|29l for the implicit scheme or in 1221 for the explicit scheme that e{n) — 0{h). We just 
have to rewrite the proof to show where constants Kg, Kf^x and Kf ^, appear precisely. For the readability of this 
article we will only give few key steps. Firstly, for the error in Y we find, for h small enough. 



sup E Yt,-Y,l ' ^ Ce 

O^k^n 



CKi 



(1 + Kl,)hl 



sup \Yt\ 

O^tsJT 



(1 + if/. J/lE 



\Zt?dt 



k=0 



\Zt 



' dt 



(A.l) 



with Zt, := iEt, Ztdt. For the error in Z we find, for h small enough, 



fc=0 













E 


sup \Yt\' 


+ 


/ \Zt\^dt 


) 




O^t^T 


Jo 





r i<t 

cj:e j 

k=0 



Zt - Zi. dt 



CK', sup E y,, - r" r .(A.2) 



OCk^n 



The study of the error X]fe=o ^ /tl''^^ I ^* ^ I '^^ done by Zhang in 1291 . Theorem 3.5 in ll26l improve 
a little bit the estimate by studying how Kg appears in the constant. Let us rewrite the proof of this theorem. 
We suppose in a first time that 6, cr, g and / are differentiable with respect to x, y and z. Then Y and Z are 
differentiable with respect to x and we obtain that 



n-l 
k=0 



\Zt-Zt \ dt 



^Ch[K^ + K]x + (1 + KjJE 



\^ZtVdt 



Thanks to classical estimates onto the solution of the BSDE solved by (VF, V Z) we have 



E 



\VZs\ ds 



^C{l + Kl+K}J{l + Kl,). 



Thus, we obtain 



fe=0 



\Zt - Zt, dt 



^Ch{l + K'^+Kl^){l+Kjj. 



(A3) 
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By standard approximation and stability results for Lipschitz BSDEs this estimate stays true when b, a, g and / 
are not differentiable. Finally, by putting together dA.ll l. (IA.2I ) and ( IA.3I ). we have 



e(n) < Che'-'^l' 



(1 + J(l + K'; + K%) + (1 + KlJE 



iZAUt 



+ (1 + KfJE 



sup \Yt\ 

0<t<T 



and the final result can be easily deduced. 



□ 
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